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Abstract of EP0990892 

Measurement data are rebinned to number of 
groups, each group has several levels parallel to 
each other and to axis of rotation. Each group 
contains one beam section (41 1-415). Rebinned 
data of each group is one-dimensional filtered in 
direction vertical to direction of levels. Spatial 
distribution of absorption is reconstruction by 
back projection of filtered data of several groups. 
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(54) Computertomographie-Verfahren mit kegelfdrmigen Strahlenbundel 



(57) Die Erfindung betrifft ein Computertomogra- 
phie-Verfahren, bei dem der Untersuchungsbereich auf 
einer Kreisbahn mit einem kegelfdrmigen Strahlenbun- 
del abgetastet wird. Der Rechenaufwand fur die Rekon- 
struktion der Absorptionsverteilung im Untersuchungs- 
bereich laBt sich dabei durch folgenden Schritte 
vereinfachen : 

a) Rebinning der MeBdaten zu einer Anzahl von 



Gruppen, wobei jede Gruppe mehrere zueinander 
und zur Rotationsachse parallele Ebenen umfaGt, 
in denen sich je ein Strahlenfacher befindet, 

b) Eindimensionale Filterung der durch das Rebin- 
ning erzeugten Daten einer jeden Gruppe in Rich- 
tung senkrecht zur Richtung der Ebenen, 

c) Rekonstruktion der raumlichen Verteilung der 
Absorption aus den gefilterten Daten verschiedener 
Gruppen. 
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Beschreibung 

[0001] Die Erfindung betrifft ein Computertomogra- 
phie-Verfahren mit den Schritten 

Erzeugen eines kegelformigen, einen Untersu- 
chungsbereich bzw. ein darin befindliches Objekt 
durchsetzenden Strahlenbundels mit einer Strah- 
lenquelle, 

Erzeugung einer kreisformigen, eine Rotation um 
eine Rotationsachse umfassenden Relativbewe- 
gung zwischen der Strahlenquelle einerseits und 
dem Untersuchungsbereich bzw. dem Objekt and- 
rerseits, 

Akquisition von MeGdaten, die von der Intensitat in 
dem Strahlenbundel jenseits des Untersuchungs- 
bereiches abhangen, mit einer Detektoreinheit 
wahrend der Relativbewegung, 
Rekonstruktion der raumlichen Verteilung der Ab- 
sorption innerhalb des Untersuchungsbereiches 
aus den akquirierten MeGdaten. 

[0002] AuGerdem bezieht sich die Erfindung auf einen 
Computertomographen zur Durchfuhrung dieses Ver- 
fahrens. 

[0003] Ein solches Verfahren ist aus dem Aufsatz von 
L.A. Feldkamp et al "Practical cone-beam algorithm", 
Journal of Optical Soc. Am.A/Vol.1 ,No 6,1984, pp. 
612-619 bekannt. Das bekannte Verfahren besteht im 
Prinzip aus den folgenden Schritten: 

a) Fur jede Strahlenquellenposition werden alle 
MeGwerte mit einem Gewichtungsfaktor multipli- 
ziert, der dem Kosinus des Winkels entspricht, den 
der Strahl langs dessen der MeGwert erfaGt wurde, 
mit dem Zentralstrahl einschlief3t. 

b) Die so gewichteten MeGwerte werden einer 
HochpaBfilterung unterzogen. 

c) Die MeGwerte werden entlang der Strahlen, 
langs derer sie gemessen worden sind in den Un- 
tersuchungsbereich ruckprojiziert. Der Beitrag ei- 
nes MeGwerts zu dem Absorptionswert eines Vo- 
xels muG dabei mit einem Faktor gewichtet werden, 
der vom Abstand des betreffenden Voxels von der 
Strahlenquellenposition abhangt. 

[0004] Da dieser letzte Schritt fur samtliche Voxel des 
zu rekonstruierenden Volumens und fur samtliche 
Strahlenquellenpositionen durchgefuhrt werden muG, 
benotigt er viel Rechenzeit. Ein weiterer Nachteil des 
bekannten Verfahrens ist, daG man damit die Absorpti- 
on nur in denjenigen Voxeln rekonstruieren kann, die 
wahrend der Untersuchung ununterbrochen von Ront- 
genstrahlung getroffen worden sind. Diese Voxel liegen 
in einem zur Rotationsachse konzentrischen diskusfor- 
migen Bereich. Erwunscht ist jedoch die Rekonstruktion 
in einem ebenen scheibenformigen Bereich. Be- 
schrankt man aber die Rekonstruktion auf den schei- 
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benformigen Bereich in dem kegelformigen Strahlen- 
bundel, der aus alien Strahlenquellenpositionen von 
Rontgenstrahlung getroffen wird, dann ergibt sich ledig- 
lich ein sehr schmaler Rekonstruktionsbereich. 

s [0005] Aufgabe der vorliegenden Erfindung ist es ein 
Com putertomographie- Verfahren der eingangs ge- 
nannten Art so auszugestalten, daG der Rechenauf- 
wand verringert wird und die Rekonstruktion der Ab- 
sorptionsverteilung in einer dickeren ebenen Scheibe 

10 des Untersuchungsbereichs moglich wird. Diese Aufga- 
be wird ausgehend von einem Verfahren der eingangs 
genannten Art erfindungsgemaG durch die folgenden 
Schritte gelost: 

15 a) Rebinning der MeGdaten zu einer Anzahl von 
Gruppen, wobei jede Gruppe mehrere zueinander 
und zur Rotationsachse parallele Ebenen umfaGt, 
in denen sich je ein Strahlenfacher befindet, 

b) Eindimensionale Filterung der durch das Rebin- 
20 ning erzeugten Daten einer jeden Gruppe in Rich- 

tung senkrecht zur Richtung der Ebenen, 

c) Rekonstruktion der raumlichen Verteilung der 
Absorption durch Ruckprojektion der gefilterten Da- 
ten von mehreren Gruppen. 

25 

[0006] Wahrend bei dem bekannten Verfahren die in 
den einzelnen Strahlenquellenpositionen gewonnenen 
MeGwerte jeweils unmittelbar einer Filterung unterzo- 
gen werden, erfolgt bei dem erfindungsgemaGen Ver- 

30 fahren zunachst ein Rebinning, bei dem in unterschied- 
lichen Strahlenquellenpositionen erfaGte Strahlenfa- 
cher gruppenweise zusammengefaGt und erst danach 
gefiltert werden. Der Vorteil des zusatzlichen Rebinning 
besteht darin, daG bei der anschlieGenden Rekonstruk- 

35 tion der Absorptionsverteilung die Beitrage der einzel- 
nen MeGwerte zu einem Voxel nicht mit einem ab- 
standsabhangigen Faktor multipliziert werden mussen. 
Dadurch vereinfacht sich das Rekonstruktionsverfahren 
ganz deutlich. Es ergibt sich auch eine verbesserte Bild- 

40 qualitat. AuGerdem ergibt sich ein gunstiges Verhaltnis 
zwischen dem bestrahlten Volumen und dem Volumen, 
innerhalb dessen die Rabsorptionsverteilung rekonstru- 
iert werden kann. 

[0007] Die Erfindung geht we iter von der Erkenntnis 
45 aus, daG sich die Absorption in jedem Voxel rekonstru- 
ieren laGt, wenn es uber einen Winkelbereich von we- 
nigstens 180° bestrahlt worden ist. Es laGt sich zeigen, 
daG alle Voxel, die diese Bedingung erfullen, in einer 
ebenen, zur Rotationsachse senkrechten Scheibe lie- 
so gen, deren Dicke groGer ist als die Abmessungen an 
den Randern des diskusformigen Bereichs, in dem die 
Absorption mit dem bekannten Verfahren rekonstruiert 
werden kann. 

[0008] Das in Anspruch 2 vorgesehene Rebinning auf 
55 einen ebenen und rechteckigen virtuellen Detektor er- 
leichtert die nachfolgende eindimensionale Filterung 
ganz wesentlich. 

[0009] In der Scheibe, in der die Absorptionsvertei- 
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lung rekonstruiert werden kann, befinden sich Voxel, die 
aus einem Winkelbereich von mehr als 180° bestrahft 
worden sind. Da fur die Rekonstruktion aber nur ein Be- 
strahlungswinkelbereich (das istder Winkelbereich, den 
die (Parallel-) Projektion der Strahlen von der Strahlen- 
quelle zu einem Voxel auf eine zur Rotationsachse 
senkrechte Ebene in dieser Ebene bedecken, bzw. den 
die Komponenten der Vektoren von der Strahlenquelle 
zu dem Voxel in der Rotationebene der Strahlenquelle 
Oberstreichen) von 180° erforderlich ist, konnen 
MeBwerte aus einigen Strahlenquellenpositionen unbe- 
rucksichtigt bleiben. Anspruch 3 beschreibt nun eine 
einfache Moglichkeit, wie man die Strahlenquellenposi- 
tionen ermitteln kann, deren MeBdaten fur die Rekon- 
struktion der Absorptionen in einzelnen Voxeln heran- 
zuziehen sind. 

[0010] Grundsatzlich kann man die Absorptionsver- 
teilung in der erwahnten Scheibe rekonstruieren, indem 
man fur alle Voxel nur so viele Strahlenquellenpositio- 
nen zur Rekonstruktion heranzieht, daB sich gerade ein 
Bestrahlungswinkelbereich von 180° ergibt. Es gibt in 
dieser Scheibe jedoch auch Voxel mit einem Bestrah- 
lungswinkelbereich von 360°. Wenn man diese Voxel 
ebenfalls nur mit einem Bestrahlungswinkelbereich von 
180° rekonstruieren wurde, wurde sich fur diese Voxel 
nicht das bestmdgliche Signal/Rauschverhaltnis erge- 
ben. Anspruch 4 beschreibt daher ein hybrides Rekon- 
struktionsverfahren, bei dem das Volumen der rekon- 
struierbaren Scheibe in zwei Teilvolumina unterteilt 
wird: in ein erstes Teilvolumen., das nur Voxel umfaBt, 
die aus alien Strahlenquellenpositionen bestrahlt wor- 
den sind und in ein zweites Teilvolumen, bei dem der 
Bestrahlungswinkelbereich kleiner ist als 360°. Bei der 
Rekonstruktion der Absorption in den Voxeln des ersten 
Teilvolumens werden samtliche MeBwerte herangezo- 
gen und bei der Rekonstruktion der Voxel in dem zwei- 
ten Teilvolumen nur die MeBwerte aus einem Bestrah- 
lungswinkelbereich von 180°. 

[0011] Anspruch 5 beschreibt einen Computertomo- 
graphen zur Durchfuhrung des erfindungsgemaBen 
Verfahrens. 

[0012] Anspruch 6 beschreibt eine Ausgestaltung des 
erfindungsgemaBen Verfahrens, die darauf abstellt, 
daB ein Teil der MeBwerte fOr die Rekonstruktion nicht 
herangezogen wird. Diese MeBwerte gehoren zu Strah- 
len, die nicht benotigt werden, wenn man nur die Ab- 
sorptionsverteilung in einer ebenen Schicht rekonstru- 
ieren will. Diese Strahlen haben in einer zur Rotations- 
achse senkrechten Ebene einen kleinen Offnungswin- 
kel und in einer die Rotationsachse enthaltenden Ebene 
einen relativ groBen Offnungswinkel. Wenn man das ke- 
gelformige Strahlenbundel durch die in Anspruch 6 an- 
gegebene Gestaltung der Kollimatoranordnung so aus- 
formt, daB diese Strahlen gar nicht erst entstehen, wird 
die Strahlenbelastung fur einen im Untersuchungsbe- 
reich befindlichen Patienten verringert. Dieser Vorteil er- 
gibt sich auch dann wenn man die Absorptionsvertei- 
lung in anderer Weise rekonstruiert, z.B. durch ein ART- 
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Verfahren. 

[0013] Die Erfindung wird nachstehend anhand der 
Zeichnungen naher erlautert. Es zeigen: 

s Fig. 1 Einen Computertomographen mit dem die Er- 
findung ausfuhrbar ist. 
Fig. 2 Ein Ablauf diagramm f Or das erf indungsgema- 

Be Rekonstruktionsverfahren. 
Fig. 3 Ein in einer Strahlenquellenposition erzeug- 
io tes kegetformiges Strahlenbundel. 

Fig. 4 Die durch das Rebinning gebildeten Strahlen- 

facher in parallelen Ebenen. 
Fig. 5 Einen Querschnitt durch diesen Strahlenfa- 
cher. 

is Fig. 6 Die geometrischen Verhaltnisse im Untersu- 
chungsbereich. 
Fig. 7 Ein Diagramm zur Auswahl von Strahlenquel- 
lenpositionen. 

20 [0014] Der in Fig. 1 dargestellte Compute rtomograph 
umfaBt eine Gantry 1 , die um eine parallel zur z-Rich- 
tung, des in Fig. 1 dargestellten Koordinatensy stems 
verlaufende Rotationsachse 14 rotieren kann. Dazu 
wird die Gantry von einem Motor 2 mit einer vorzugs- 

2s weise konstanten, aber einstellbaren Winkelgeschwin- 
digkeit angetrieben. An der Gantry ist eine Strahlenquel- 
le S, beispielsweise ein Rontgenstrahler, befestigt. Die- 
ser ist mit einer Kollimatoranordnung 3 versehen, die 
aus der von der Strahlenquelle S erzeugten Strahlung 

30 ein kegetformiges Strahlenbundel 4 ausblendet, d.h. ein 
Strahlenbundel, das sowohl in z-Richtung als auch in 
einer dazu senkrechten Richtung (d.h. in einer zur Ro- 
tationsachse senkrechten Ebene) eine von Null ver- 
schiedene, endliche Ausdehnung hat. 

35 [001 5] Das Strahlenbundel 4 durchdringt einen Unter- 
suchungsbereich 1 3, in dem sich ein Patient auf einem 
Patientenlagerungstisch (beides nicht naher darge- 
stellt) befinden kann. Der Untersuchungsbereich 13 hat 
die Form eines Zylinders, der im folgenden als Objekt- 

40 zylinder 13 bezeichnet wird. Nach dem Durchsetzen 
des Objektzylinders 1 3 trifft das Rontgenstrahlenbundel 
4 auf eine an der Gantry 1 befestigte zweidimensionale 
Detektoreinheit 16, die eine Anzahl von Detektorzeilen 
mit jeweils einer Vielzahl von Detektorelementen um- 

45 faBt. Die Detektorzeilen befinden sich in zur Rotations- 
achse senkrechten Ebenen auf einem Kreisbogen um 
die Strahlenquelle S; sie konnen aber auch einen Kreis- 
bogen um die Rotationsachse beschreiben oder gerad- 
linig sein. Jedes Detektorelement liefert in jeder Strah- 

50 lenquellenposition einen MeBwert fur einen Strahl aus 
dem StrahlenbOndel 4. 

[001 6] Der mit a max bezeichnete Offnungswinkel des 
Strahlenbundels 4 (als Offnungswinkel ist der Winkel 
definiert, den ein Strahl den einer zur Rotationsachse 
55 14 senkrechten Ebene am Rande des Strahlenbundels 
4 liegt, mit einer durch die Strahlenquelle S und die Ro- 
tationsachse 14 definierten Zentralstrahlebene ein- 
schlieBt) bestimmt dabei den Durchmesser des Objekt- 
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zylinders 1 3, innerhalb dessen sich das zu untersuchen- 
de Objekt bei der Akquisition der MeGwerte befindet. 
Der Untersuchungsbereich 13 - bzw. das Objekt oder 
der Patientenlagerungstisch - kann mittels eines Motors 
5 parallel zur Rotationsachse 14 bzw. zur z-Achse ver- 
schoben werden. Die Geschwindigkeit dieses Vor- 
schubs in z-Richtung ist vorzugsweise konstant und ein- 
stellbar. Wenn die Motoren 5 und 2 gleichzeitig laufen, 
ergibt sich eine helixformige Abtastbewegung der 
Strahlenquelle S und der Detektoreinheit. Wenn hinge- 
gen der Motor 5 fur den Vorschub in z-Richtung stillsteht 
und der Motor 2 die Gantry einzeln rotieren laGt, ergibt 
sich eine kreisformige Abtastbewegung der Strahlen- 
quelle S und der Detektoreinheit relativ zum Untersu- 
chungsbereich 1 3. Im folgenden soil nur diese kreisfor- 
mige Abtastbewegung betrachtet werden. 
[0017] Die von der Detektoreinheit akquirierten 
MeGdaten werden einem Bildverarbeitungsrechner 10 
zugefuhrt, der daraus die Absorptionsverteilung in ei- 
nem Teil des Untersuchungsbereichs 13 rekonstruiert 
und z.B. auf einem Monitor 11 wiedergibt. Die beiden 
Motoren 2 und 5, der Bildverarbeitungsrechner 10, die 
Strahlenquelle S und der Transfer der MeGdaten von 
der Detektoreinheit 16 zum Bildverarbeitungsrechner 
10 werden von einer geeigneten Kontrolleinheit gesteu- 
ert. 

[001 8] Fig. 2 zeigt den Ablaut eines MeG- und Rekon- 
struktionsverfahrens, das mit dem Computertomogra- 
phen nach Fig. 1 durchgefuhrt werden kann. 
[0019] Nach der Initialisierung im Block 101 rotiertdie 
Gantry mit einer konstanten Winkelgeschwindigkeit. Im 
Schritt 102 wird die Strahlung der Strahlenquelle S ein- 
geschaltet, und die dabei von den Detektorelementen 
der Detektoreinheit 16 akquirierten MeGwerte werden 
in einem Speicher des Bildverarbeitungsrechners 10 
gespeichert 

[0020] Fig. 3 zeigt die kreisformige Bahn 17, auf der 
sich die Strahlenquelle S und die Detektoreinheit 16 um 
die Rotationsachse 14 bewegen. Fur eine bestimmte 
Strahlenquellenposition S 0 ist das Strahlenbundel 4 
dargestellt Man kann sich dieses kegelformige Strah- 
lenbundel 4 aus einer Vielzahl von ebenen Strahlenfa- 
chern zusammensetzen, die sich - wie die Strahlenfa- 
cher 401 ...403 in zur Rotationsachse parallelen Ebe- 
nen befinden. Obwohl man die Strahlen des kegelfor- 
migen Strahlenbundels 4 auch auf andere Weise zu- 
sammenfassen kann, soil der Begriff "Strahlenfacher" 
im folgenden nur fur solche Strahlen verwendet werden, 
die in einer gemeinsamen, zur Rotationsachse 14 par- 
allelen Ebene liegen. Diese Strahlenfacher gehen von 
einer Strahlenquellenposition aus und werden von je ei- 
ner zur Rotationsrichtung parallelen Spalte von Detek- 
torelementen auf der Detektoreinheit 16 erfaGt. Fig. 3 
deutet an, daG auch in anderen Positionen der Strah- 
lenquelle (z.B. S.-,, St Oder S 2 ) das emittierte kegelfor- 
mige Strahlenbundel gemessen wird. 
[0021] Im Schritt 1 03 wird der Durchmesser des dar- 
zustellenden Bereichs FOV (field ofview) festgelegt. Der 



Durchmesser kann dem durch definierten Durch- 
messer des Untersuchungsbereichs 1 3 entsprechen; er 
kann jedoch auch kleiner sein. 

[0022] In den Schritten 104-107 erfolgt ein Rebinning. 

5 Die Daten werden dabei so umsortiert und uminterpo- 
liert, als waren Sie mit einer anderen Strahlenquelle (ei- 
ner kreisbogenformigen Strahlenquelle, die zueinander 
parallele Strahlenfacher emittiert) und mit einem ande- 
ren Detektor (einem ebenen, rechteckigen und die Ro- 

10 tationsachse enthaltenden ,,virtuellen B Detektor) ge- 
messen worden. 

[0023] Im Schritt 104 werden dazu zunachst die 
Strahlenfacher von verschiedenen Strahlenquellenpo- 
sitionen zu jeweils einer Gruppe zusammengefaGt, die 

15 in zueinander parallelen Ebenen liegen. Die zu einer 
Gruppe gehorenden Strahlenfacher erfullen daher die 
Bedingung, daG die Summe der Winkel a und (5 fur alle 
Strahlenfacher dieser Gruppe den gleichen Wert hat. 
Dabei ist a der Winkel, den die Ebene des Strahlenfa- 

20 chers mit einer durch die Strahlenquellenposition und 
die Rotationsachse definierten Zentralstrahlebene ein- 
schlieGt und der durch die Lage der Spalte von Detek- 
torelementen gegeben ist, die den betreffenden Strah- 
lenfacher erfaGt hat. p ist ein die Strahlenquellenpositi- 

25 on (z.B. S 0 ) auf dem Kreis 17 kennzeichnender Winkel. 
Wenn die Strahlenfacher fur eine Strahlenquellenposi- 
tion diese Bedingung nicht exakt erfullen, muG fur diese 
Strahlenquellenposition ein entsprechender Strahlenfa- 
cher durch Interpolation aus den Strahlen benachbarter 

30 Strahlenfacher ermittelt werden. 

[0024] Fig. 4 zeigt eine auf diese Weise gebildete 
Gruppe von Strahlenfachern. Von jeder der einander 
benachbarten Strahlenquellenpositionen S_ 2 ....S 0 ,...S 2 
gehort jeweils ein Strahlenfacher zu einer Gruppe. Jede 

35 Gruppe kann durch den Winkel p der Strahlenquellen- 
position (z.B. S 0 ) gekennzeichnet werden, deren zu der 
Gruppe gehorender Strahlenfacher die Rotationsachse 
14 durchsetzt. (In der Regel ist dies die mittlere Strah- 
lenquellenposition - im Beispiel nach Fig. 4 also die 

40 Strahlenquellenposition S 0 ). - Es konnen sich dann so 
viele Gruppen von Strahlenfachern ergeben , wie es 
Strahlenquellenpositionen gibt, aber auch mehr oder 
weniger. 

[0025] Die auf diese Weise bestimmten Strahlenfa- 
45 cher - darunter die in Fig. 4 dargestellten Strahlenfacher 
411. ..41 5 - definieren ein Strahlenbundel 410, das eine 
zeltartige Form hat und sich aus Strahlenfachern zu- 
sammensetzt, die in zueinander und zur Rotationsach- 
se parallelen Ebenen liegen. In Fig. 4 ist auGerdem der 
50 Schnittbereich 420 dargestellt, der sich ergibt, wenn das 
StrahlenbOndel 410 von einer die Rotationsachse 14 
enthaltenden und zu den Ebenen der Strahlenfacher 
411 ...415 senkrechten Ebene geschnitten wird. Wie ins- 
besondere Fig. 5 zeigt, die diesen Schnittbereich dar- 
55 stellt, sind die oberen und unteren Rander gewdlbt. Die- 
se Wolbung ist darauf zuruckzufuhren, daG die Strah- 
lenquellen positionen in der Mitte (z.B. S 0 ) weiter von 
der Schnittebene entf ernt sind als die am Rande S 2 oder 
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S_ 2 und daB die Facher alle den gleichen Offnungswin- 
kel haben, weil die Detektorzeilen einen Kreisbogen urn 
die Strahlenquelle S beschreiben. Bei einer anderen 
Geometrie der Detektorzeilen wurde sich die Form der 
Schnittebene 420 andern. Bei einer ebenen Detektor- 
einheit (z.B. mit gradlinigen Detektorzeilen) ware die 
Wolbung noch ausgepragter, weil die am Rande liegen- 
den Strahlenfacher , wie z.B. 411 und 415, dabei einen 
geringeren Offnungswinkel hatten. 
[0026] Fur jede Gruppe von Strahlenfachem wird da- 
her in dem ebenen Schnittbereich 420 ein rechteckiger 
virtueller Detektor definiert (Schritt 105), dessen oberer 
und unterer Rand 161 bzw. 162durch die Abmessungen 
der auBeren Strahlenfacher 41 1 bzw. 41 5 in den ebenen 
Schnittbereich gegeben sind. Wenn der Durchmesser 
des zu rekonstruierenden FOV (field of view) kleiner ge- 
wahlt worden ist als der Durchmesser des Untersu- 
chungsbereichs 13, konnen von jeder Gruppe jeweils 
die auBeren Strahlenfacher (z.B. 411, 415) entfallen. 
Der Abstand der oberen und unteren Rander 161 und 
162 kann in diesem Fall groBer gewahlt werden als bei 
dem in Fig. 5 dargestellten Beispiel. 
[0027] In Fig. 5 sind auBerdem - durch runde Punkte 
markiert - die DurchstoBpunkte einiger in den Strahlen- 
fachem 411—415 enthaltener Strahlen durch diesen 
virtuellen Detektor dargesteilt. SchlieBlich sind durch 
Kreuze die Stutzpunkte eines regelmaBigen kartesi- 
schen Gitters angedeutet. Die DurchstoBpunkte und die 
Stutzpunkte fallen durchweg nicht zusammen. Man er- 
kennt einerseits, daB die Strahlenfacher auBen enger 
beieinanderliegen als innen und daB die 
DurchstoBpunkte eines Strahlenfachers innen weiter 
voneinander entfernt sind als auBen. Aus den MeBwer- 
ten fur die DurchstoBpunkte mussen daher in den bei- 
den folgenden Schritten 106 und 107 die MeBwerte an 
den aquidistanten Stutzstellen innerhalb des virtuellen 
Detektors 160 ermittelt werden. 

[0028] Im Schritt 1 06 erfolgt dabei zunachst eine ver- 
tikale Interpolation derart, daB die Stutzstellen in verti- 
kaler Richtung fur alle Strahlenfacher voneinander den 
gleichen Abstand haben, wie die DurchstoBpunkte bzw. 
Stutzstellen am Rande des virtuellen Detektors 160. Im 
Schritt 107 erfolgt eine Interpolation in horizontaler 
Richtung, so daB man interpolierte Werte f Or Stutzstel- 
len erhalt, die innerhalb des virtuellen Detektors 160 in 
horizontaler Richtung voneinander den gleichen Ab- 
stand haben. Die Interpolationsschritte 105 und 106 
konnen auch miteinander vertauscht und ggf. sogar mit- 
einander kombiniert werden. 

[0029] Das Ergebnis dieser Interpolation in horizon- 
taler und vertikaler Richtung ist, daB innerhalb des vir- 
tuellen Detektors 1 60 die Intensitat der Strahlung an den 
aquidistanten Stutzpunkten eines regelmaBigen karte- 
sischen Gitters vorliegt. Die Stutzstellen definieren ih- 
rerseits wiederum neue Strahlenfacher, die in parallelen 
Ebenen mit jeweils dem gleichen Abstand voneinander 
liegen. Nur diese neuen Strahlenfacher (die teilweise 
mit den originalen Strahlenfachem identisch sein kon- 



nen) werden fur die weitere Rekonstruktion benutzt 
[0030] Ein Teil der Strahlen, deren DurchstoBpunkte 
auBerhalb des virtuellen Fensters liegen, werden in den 
Schritten 106 und 107 nicht benutzt. Es ist daher von 

s Vorteil, den Kollimator 3 von vornherein so zu gestalten, 
daB das kegelfdrmige Strahlenbundel von vornherein 
nicht enthalt. Dadurch wurde die Strahlenbelastung fur 
den Patienten verringert. Anstelle von geradlinigen 
senkrecht zur Rotationsachse verlaufenden Kanten 

10 muBte die Kollimatoranordnung zu diesem Zweck nach 
innen gewolbte Kanten aufweisen, so daB die Strahlen- 
facher, die die Rotationsachse schneiden oder davon 
einen geringen Abstand haben, einen geringeren Off- 
nungswinkel hatten, als Strahlenfacher am auBeren 

is Rand des Strahlenbundels. 

Nach den Schritten 103 bis 107 ist somitfur jede Gruppe 
von Strahlenfachem die Strahlenintensitat an den regel- 
maBigen Stutzstellen des zu dieser Gruppe gehorenden 
virtuellen Detektors 160 bestimmt. Dies erleichtert die 

20 erforderliche HochpaBfilterung fundamental, weil ledig- 
lich eine eindimensionale Filterung der durch das Re- 
binning auf den virtuellen Detektor 1 60 erzeugten Daten 
notwendig ist. Im Schritt 108 wird daher auf diese Daten 
eine eindimensionale Filterung mit einem rampenformig 

25 mit der Frequenz ansteigenden Ubertragungsfaktor an- 
gewandt. Dazu werden jeweils lediglich die Werte von 
in horizontaler Richtung aufeinander folgenden Stutz- 
stellen herangezogen. Diese Filterung wird fur alle 
Gruppen von (neuen) Strahlenfachem ausgefuhrt. 

30 [0031] Die nach dem Rebinning und der Filterung fur 
die durch die Stutzstellen im virtuellen Fenster 160 de- 
finierten Strahlen ermittelten Daten werden anschlie- 
Bend zur Rekonstruktion der Absorptionsverteilung im 
Untersuchungsbereich durch eine Ruckprojektion her- 

^5 angezogen. Durch das vorangegangenen Schritte 
103-108 ist es dabei aber - im Gegensatz zu dem ein- 
gangs erwahnten bekannten Verfahren - nicht erforder- 
lich, die Beitrage einzelner Strahlen zu dem Absorpti- 
onswert eines Voxels mit einem Faktor zu gewichten 

40 werden, der vom Abstand des betreffenden Voxels von 
der Strahlenquellenposition abhangt. 
[0032] Im Prinzip wurde es fur die Ruckprojektion ge- 
nugen, wenn man fur jedes Voxel Daten aus Gruppen 
heranzieht, deren Strahlenfacher miteinander einen 

45 Winkelbereich von insgesamt nur 180° definieren. Es 
gibt jedoch im Untersuchungsbereich eine Vielzahl von 
Voxeln, die in samtlichen Strahlenquellenpositionen von 
Strahlung getroffen worden sind und es ist empfehlens- 
wert fur diese Voxel samtliche verfugbaren MeBwerte 

50 heranzuziehen, um ein moglichst gutes Signal/Rausch- 
verhattnis zu erreichen. Deshalb wird im folgenden da- 
nach differenziert, ob die Voxel aus alien Strahlenquel- 
lenpositionen bestrahlt worden sind oder nicht. 
[0033] Nach der Vorgabe eines Voxels im Untersu- 

55 chungsbereich, dessen Absorption rekonstruiert wer- 
den soli (Schritt 109), wird im Schritt 110 zunachst ge- 
pruft, ob das betreffende Voxel aus alien Strahlenquel- 
lenpositionen bestrahlt worden ist. 
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[0034] Dazu wird zunachst auf Fig. 6 hingewiesen, die 
die durch einen Punkt angedeutete Strahlenquelle S, 
die durch eine Linie symbolisierte Detektoreinheit 16 
und das Strahlenbundel 4 in einer ersten Position be- 
zOglich der Rotationsachse 14 zeigt. AuBerdem zeigt 
Fig. 6 in einer zweiten, gegenuber der ersten um 180° 
versetzten Position die Strahlenquelle S\ die Detektor- 
einheit 16' und das Strahlenbundel 4'. AuBerdem ist in 
Fig. 6 der Untersuchungsbereich 13 angedeutet. Die 
Strahlen an den in Richtung der Rotationsachse 14 ge- 
geneinander versetzten Randern des Strahlenbundels 
4 sind gestrichelt angedeutet. Mit ausgezogenen Linien 
sind hingegen die Strahlen 40 bzw. 40' angedeutet, die 
nach dem Weglassen der Strahlen am Rande (Schritt 
107) das Strahlenbundel begrenzen. 
[0035] Es laBt sich zeigen, daB alle Voxel innerhalb 
einer ebenen kreisformigen Scheibe V 0 aus einem Win- 
kelbereich von zumindest 180° bestrahlt worden sind. 
Die Seitenflachen dieser Scheibe sind durch die Schnitt- 
punkte der Randstrahlen 40 bzw. 40' mit der Rotations- 
achse 14 definiert. Alle Voxel, die innerhalb des durch 
die Randstrahlen 40 bzw. 40' definierten, diskusformi- 
gen Teilvolumens V t liegen, sind dabei aus samtlichen 
Strahlenquellenpositionen bestrahlt worden, wahrend 
die Voxel, die in dem durch die Differenz von V 0 und V n 
gegebenen Teilvolumen V 2 liegen, die Strahlenquelle 
bei der Abtastbewegung aus einem Winkel von minde- 
stens 180° aber weniger als 360° "gesehen" haben. 
[0036] Ergibt nun die Abf rage 1 1 0, daB das im Schritt 
109 ausgewahlte Voxel im Teilvolumen V-, liegt, daB es 
also aus alien Strahlenquellenpositionen bestrahlt wor- 
den ist, dann erfolgt im Schritt 111 eine ROckprojektion, 
wobei von alien Gruppen jeweils derjenige Strahi in die- 
ser Gruppe herangezogen wird, der durch das betref- 
fende Voxel verlauft. Passiert kein Strahi exakt die Mitte 
dieses Voxels, dann muB der zugehorige Wert durch In- 
terpolation der MeBwerte von benachbarten Strahlen 
ermittelt werden. 

[0037] Ist im Schritt 110 die Abf rage negativ, dann er- 
folgt im Schritt 112 eine weitere Abfrage, ob das ausge- 
wahlte Voxel im Teilvolumen V 2 liegt (ob es also die 
Strahlenquelle unter einem Winkel von wenigstens 180° 
"gesehen" hat) ist dies nicht der Fall, dann verzweigt das 
Programm zu einer weiteren Abfrage 11 3, in der festge- 
stellt wird, ob die Absorption schon in samtlichen Voxeln 
des FOV rekonstruiert worden ist. Ist dies nicht der Fall, 
dann werden die Schritte 109 ff erneut durchlaufen. Bei 
einer positiven Ant wort, d.h. wenn sich das im Schritt 
109 ausgewahlte Voxel in dem Volumen V 2 befindet, 
werden im Schritt 118 die Gruppen ausgewahlt, die zur 
Rekonstruktion der Absorption in dem Voxel herange- 
zogen werden. Dazu wird auf Fig. 7 verwiesen, die die 
kreisformige Abtastbahn 17 zeigt und auf ihr - durch 
Punkte angedeutet - die verschiedenen Strahlenquel- 
lenpositionen. Die in Fig. 4 mit S_ 2 ...S 0 ...S 2 bezeichne- 
ten Strahlenquellenpositionen sind in Fig. 7 durch die 
zugehorigen Winkel p_ 2 ..,p 0 ...p 2 auf dem Kreis 17 ge- 
kennzeichnet. Wie in Verbindung mit den Schritten 103 



bis 108 erlautert, ist jeder Strahlenquellenposition eine 
Gruppe von durch das Rebinning entstandenen gefilter- 
ten MeBdaten zugeordnet 

[0038] AuBerdem ist in Fig. 7 ein Voxel P 2 dargestellt, 

5 das - wie in Fig. 6 angedeutet - in dem Teilvolumen V 2 
liegt. Mit E 2 ist dabei die zur Zeichenebene der Fig. 7 
senkrechte Ebene bezeichnet, die durch das Voxel P 2 
und die Rotationsachse 14 definiert wird. SchlieBlich ist 
mit E eine Ebene bezeichnet, die die Ebene E 2 in der 

io Rotationsachse schneidet. 

[0039] Es laBt sich zeigen, daB das Voxel P 2 aus alien 
Strahlenquellenpositionen jenseits der Ebene E be- 
strahlt worden ist. Hingegen hat das Voxel P 2 allenfalls 
von einigen - nicht aber von alien - diesseits der Ebene 

15 e (in bezug auf das Voxel) befindlichen Strahlenquel- 
lenpositionen Strahlung empfangen. 
[0040] Infolgedessen werden im Schritt 1 1 9 alle durch 
das Rebinning und das Filtern entstandenen Gruppen 
von Daten zur Rekonstruktion herangezogen, die den 

20 jenseits der Ebene E befindlichen Strahlenquellenposi- 
tionen (z.B. p 2 ...0. 2 ) zugeordnet sind und die miteinan- 
der einen Bestrahlungswinkelbereich von genau 180° 
bilden. Dieser Bestrahlungswinkelbereich ist notwendig 
und hinreichend, um die Absorption in dem betreffenden 

25 Voxel zu rekonstruieren, auch wenn die mit diesen Da- 
ten durchgefOhrte ROckprojektion ein schlechteres Si- 
gnal/Rauschverhaltnis ergibt als die im Schritt 111 
durchgefOhrte ROckprojektion. 

[0041] Sobald im Schritt 113 festgestellt ist, daB alle 
30 Voxel in dem ausgewahlten FOV rekonstruiert worden 
sind, ist das Verfahren damit beendet (Schritt 114). Die 
Absorptions werte in dem FOV konnen dann - ggf . schei- 
benweise - auf einem geeigneten Display wiedergege- 
ben werden. 

35 

Patentanspruche 

1 . Computertomographie-Verfahren mit den Schritten 

40 

Erzeugen eines kegelformigen, einen Untersu- 
chungsbereich bzw. ein darin befindliches Ob- 
jekt durchsetzenden Strahlenbundels mit einer 
Strahlenquelle, 

45 - Erzeugung einer kreisformigen, eine Rotation 
um eine Rotationsachse umfassenden Relativ- 
bewegung zwischen der Strahlenquelle einer- 
seits und dem Untersuchungsbereich bzw. 
dem Objekt andrerseits, 

50 - Akquisition von MeBdaten, die von der Intensi- 
tat in dem Strahlenbundel jenseits des Unter- 
suchungsbereiches abhangen, mit einer De- 
tektoreinheit wahrend der Relativbewegung, 
Rekonstruktion der raumlichen Verteilung der 

55 Absorption innerhalb des Untersuchungsberei- 

ches aus den akquirierten MeBdaten 

gekennzeichnet durch folgende Schritte: 
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a) Rebinning der MeGdaten zu einer Anzahl 
von Gruppen, wobei jede Gruppe mehrere zu- 
einander und zur Rotationsachse parallels 
Ebenen umfaGt, in denen sich je ein Strahlen- 
facher befindet, 5 

b) Eindimensionale Filterung der durch das Re- 
binning erzeugten Daten einer jeden Gruppe in 
Richtung senkrecht zur Richtung der Ebenen, 

c) Rekonstruktion der raumlichen Verteilung 
der Absorption durch Ruckprojektion der gefil- 10 
terten Daten von mehreren Gruppen. 

2. Computertomographie-Verfahren nach Anspruch 
1, 

dadurch gekennzeichnet, daG das Rebinning auf je is 
einen zu den Ebenen einer jeden Gruppe senkrech- 
ten, ebenen und rechteckigen virtuellen Detektor 
erfolgt, der die Rotationsachse enthalt. 

6. 

3. Computertomographie-Verfahren nach Anspruch 20 
1, 

gekennzeichnet durch die Rekonstruktion der Ab- 
sorption der Voxel mit einem Bestrahlungswinkel- 
bereich von mindestens 180° unter Heranziehung 
der gefilterten MeGdaten von ausschlieGlich denje- 25 
nigen Strahlenquellenpositionen, die von dem Vo- 
xel durch einedie Rotationsachse enthaltende Ebe- 
ne getrennt sind, welche zu einer durch das betref- 
fende Voxel und die Rotationsachse definierten 
Ebene senkrecht steht. 30 

4. Computertomographie-Verfahren nach Anspruch 
1, 

gekennzeichnet durch folgende Schritte: 

35 

a) Zuordnung von Voxeln mit einem Best rah - 
lungswinkelbereich von 360° zu einem ersten 
Teilvolumen, 

b) Zuordnung von Voxeln mit einem Bestrah- 
lungswinkelbereich von mindestens 180° aber 40 
weniger als 360° zu einem zweiten Teilvolu- 
men, 

c) Rekonstruktion der Absorption der Voxel in 
dem ersten Teilvolumen unter Heranziehung 
von MeGdaten von samtiiche Strahlenquellen- 45 
positionen, 

d) Rekonstruktion der Absorption der Voxel in 
dem zweiten Teilvolumen unter Heranziehung 
der gefilterten MeGdaten von gerade so vielen 
Strahlenquellenpositionen, , daG sich ein Be- 50 
st rah lungswinkelbereich von 180° ergibt.. 

5. Computertomograph zur Durchfuhrung des Verfah- 
rens nach Anspruchs 1 mit einer Strahlenquelle und 
einer damit gekoppelten Detektoreinheit sowie mit 55 
einer Antriebsanodnung, um ein im Untersuchungs- 
bereich enthaltenes Objekt und die relativ zueinan- 
der um eine Rotationsachse rotieren zu lassen und 



mit einer Rekonstruktionseinheit zur Rekonstrukti- 
on der raumlichen Verteilung der Absorption inner- 
halb des Untersuchungsbereiches aus den von der 
Detektoreinheit akquirierten MeGdaten, 
gekennzeichnet durch 

a) Mittel zum Rebinning der MeGdaten zu meh- 
reren Gruppen, wobei jede Gruppe mehrere 
zur Rotationsachse parallele Ebenen umfaGt, 
in denen sich je ein Strahlenfacher befindet, 

b) Mittel zur eindimensionalen Filterung der 
durch das Rebinning erzeugten Daten einer je- 
den Gruppe in Richtung senkrecht zu zur Rich- 
tung den Ebenen 

c) Mittel zur Rekonstruktion der raumlichen 
Verteilung der Absorption aus den gefilterten 
Daten verschiedener Gruppen 

Computertomograph, insbesondere nach An- 
spruch 5, 

dadurch gekennzeichnet, daG die in Richtung der 
Rotationsachse gegeneinander versetzten Rander 
der Kollimatoranordnung so geformt sind, daG das 
kegelformige Strahlenbundel in seiner Mitte eine 
geringere Offnung hat als an seinen Randern. 
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FIG.3 




FIG. A- 
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FIG.7 
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